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ABSTRACT 

We have modeled nova light curves exceeding the Eddington luminosity. It has been suggested that a 
porous structure develops in nova envelopes during the super Eddington phase and the effective opacity 
is much reduced for such a porous atmosphere. Based on this reduced opacity model, we have calculated 
envelope structures and light curves of novae. The optically thick wind model is used to simulate 
nova winds. We find that the photospheric luminosity and the wind mass- loss rate increase inversely 
proportional to the reducing factor of opacities, but the wind velocity hardly changes. We also reproduce 
the optical light curve of V1974 Cygni (Nova Cygni 1992) in the super-Eddington phase, which lasts 13 
days from the optical peak 1.7 mag above the Eddington luminosity. 

Subject headings: stars: individual (V1974 Cyg) — stars: interior — stars: mass loss — novae — stars: 
white dwarfs 



1. INTRODUCTION 

The super-Eddington luminosity is one of the long 
standing problems in nova theory. The peak luminosity 
of classical novae often exceeds the Eddington limit by a 
factor of several (Delia Valle & Livio 1995, and refer- 

. ences therein). Super-Eddington phases last several days 
or more, longer than the dynamical time scale of white 

' dwarf (WD) envelopes. Many theoretical works have at- 
tempted to reproduce this phenomena, but not succeeded 
yet (e.g., Prialnik et al 1978; Sparks et al 1978; Nariai et 
al. 1980; Starrfield et al. 1985; Prialnik et al 1986; Kutter 
& Sparks 1989; Prialnik & Kovetz 1992; Kovetz 1998). 

The Eddington luminosity is understood as an upper 
limit of the luminosity of stars in hydrostatic equilibrium. 
Its classical definition is 
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where c is the light speed, G the gravitational constant, 
M the mass of the WD, and K e i is the opacity by electron 
scattering. If the diffusive luminosity exceeds this limit, 
the stellar envelope cannot be in hydrostatic balance and 
a part of the envelope is ejected. 

During nova outbursts nuclear burning produces energy 
much faster than this limit. The envelope is accelerated 
deep inside the photosphere and a part of it is ejected as 
a wind. Once the wind occurs, the diffusive luminosity is 
consumed to drive the wind. As a result, the photospheric 
luminosity decreases below the Eddington limit defined by 
equation (1) (Kato 1983, 1985). 

Recently, Shaviv (2001b, 2002) presented a new idea 
of clumpy atmospheres to explain the super-Eddington lu- 
minosity of novae. The nova envelope becomes unstable 
against dumpiness shortly after the ignition when the lu- 



minosity exceeds a critical fraction of the Eddington limit 
(Shaviv 2001a). Such a clumpy structure reduces effec- 
tive opacity, and correspondingly, increases the effective 
Eddington luminosity. Therefore, the luminosity could be 
larger than the classical Eddington limit, even though it 
does still not exceed the effective Eddington luminosity. 

Shaviv (2001b) suggested a model of nova envelope with 
the super-Eddington luminosity consisting of four parts: 
(1) convective region: a bottom region of the envelope in 
which the diffusive luminosity is sub-Eddington and addi- 
tional energy is carried by convection, (2) a porous atmo- 
sphere: the effective Eddington luminosity is larger than 
the classical Eddington limit, (3) an optically thick wind 
region: the effective Eddington limit tends to the classical 
value, and (4) the photosphere and above. 

Based on Shaviv's picture, we have assumed reduced 
opacities to model the super-Eddington phase of V1974 
Cyg (Nova Cygni 1992). V1974 Cyg is a well observed 
classical nova so that various multiwavclcngth observations 
are available such as optical (Iijima and Esenoglu 2003), 
supersoft X-ray (Krautter 2002), and radio (Eyres at al. 
2005). Chochol et al. (1997) summarized observational 
estimates of optical maximum magnitude, ranging from 
—7.3 to —8.3 mag with an average magnitude of —7.78. 
These values indicate that the peak luminosity exceeded 
the Eddington limit by more than a magnitude and the 
duration of super-Eddington phase lasts several days or 
more. 

In §2, we briefly describe our numerical method. Physi- 
cal properties of the envelope with reduced effective opac- 
ities are shown in §3. Our light curve model of V1974 Cyg 
is given in §4. Discussion follows in §5. 

2. ENVELOPE MODEL WITH REDUCED OPACITY 
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We have calculated structures of envelopes on mass- 
accreting WDs by solving the equations of motion, mass 
continuity, energy generation, and energy transfer by diffu- 
sion. The computational method and boundary conditions 
are the same as those in Kato & Hachisu (1994) except 
the opacity. We use an arbitrarily reduced opacity 

K eS = K/s, (2) 

where k is OPAL opacity (Iglesias & Rogers 1996) and s 
is a opacity reduction factor that represents reduced ratio 
of opacity due to dumpiness of the envelope. The effective 
Eddington luminosity now becomes 
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When s is greater than unity, the luminosity can be larger 
than the classical Eddington limit (1). Note that the Ed- 
dington luminosity (3) is a local variable because OPAL 
opacity is a function of local variables. 

As a first step, we simply assume that the opacity re- 
duction factor s is spatially constant. Figure 1 shows nu- 
merical results for three envelopes of s = 1, 3, and 10 
on a 1.0 Mq WD with the Chandrasckhar radius, i.e., 
logi?wD (cm) = 8.733. The chemical composition of 
the envelope is assumed to be uniform, i.e., X = 0.35, 
Y = 0.33, C + 0=0.3, and Z = 0.02, where Z includes 
carbon and oxygen by solar composition ratio for heavy 
elements. 

In Figure 1, the effective Eddington luminosity (3) 
is plotted by dashed lines, which sharply decreases at 
logr (cm) ~ 11.1 corresponding to the iron peak in OPAL 
opacity at logT (K) ~ 5.2. The wind is accelerated in 
this region and reaches a terminal velocity deep inside the 
photosphere. The diffusive luminosity (L r ) decreases out- 
ward because the energy flux is consumed to push matter 
up against the gravity. These features are qualitatively 
the same as those in the three nova envelopes with s = 1, 
3, and 10. 

Figure 2 shows the photospheric velocity (v p h), the wind 
mass loss rate (M), and the photospheric luminosity (T p h) 
for three evolutionary sequences of s = 1,3, and 10. The 
s = 1 sequence is already reported in Kato & Hachisu 
(1994). In each evolutionary sequence, the envelope mass 
is large for smaller photospheric temperature (T p h). The 
figure also shows that T p h and M increase almost pro- 
portionally to s, whereas the wind velocity (v p h) hardly 
changes but even slightly decreases. 

Theoretical light curves are calculated from these se- 
quences. After the onset of a nova outburst, the enve- 
lope expands to a giant size and the luminosity reaches its 
peak. After that, the envelope mass gradually decreases 
owing mainly to the wind mass loss. During the nova 
decay phase, the bolometric luminosity is almost constant 
whereas the photospheric temperature increases with time. 
The main emitting wavelength region moves from optical 
to supersoft X-ray through ultra-violet (UV). Therefore, 
we obtain decreasing visual magnitudes (Kato & Hachisu 
1994). 

Figure 3 shows visual light curves for the opacity reduc- 
tion factor s = 1, 3, and 10. The visual magnitude decays 
quickly for a larger s because an envelope for a larger s 



has a heavier wind mass loss. The peak luminosity of 
each light curve is shown by arrows. When the opacity 
reduction factor s is larger than unity, the peak luminos- 
ity exceeds the classical Eddington limit, which is roughly 
corresponding to the Eddington luminosity for s = 1. 

3. LIGHT CURVE OF V1974 CYG 

Recently, Hachisu & Kato (2005) presented a light curve 
model of V1974 Cyg that reproduced well the observed X- 
ray, UV, and optical light curves except for a very early 
phase of the super-Eddington luminosity. Here, we focus 
on this early phase (m v > 6) and reproduce the super- 
Eddington luminosity based on the reduced opacity model. 

We adopt various WD model parameters after their 
best fit model, i.e., the mass of 1.05 M Q , radius of 
log (R/Rq) = —2.145, and chemical composition of 
X=0A6, CNO = 0.15, iVe = 0.05, and Z = 0.02 by mass. 
These parameters are determined from the X-ray turn-off 
time, epoch at the peak of UV 1455 A flux, and epoch at 
the wind termination. All of these epochs are in the post 
super-Eddington phase. 

Our simple model with a constant s such as in Figure 3 
does not reproduce the observed light curve of V1974 Cyg. 
Therefore, we assumed that s is a decreasing function of 
time. Here, the decreasing rate of s is determined from the 
wind mass loss rate and the envelope mass of solutions we 
have chosen. After many trials, we have found that we 
cannot obtain a light curve as steep as that of V1974 Cyg. 
Finally, we further assume that s is a function both of tem- 
perature and time. We define s as unity in the outer part 
of the envelopes (logT < 4.7), but a certain constant value 
(s > 1) in the inner region (logT > 5.0), and changes lin- 
early between them. This assumption well represents the 
nova envelope model by Shaviv (2001b) outlined in §1. 
After many trials, we choose s — 5.5 at the optical peak 
(JD 2,448,676) and gradually decreases it to 1.0 with time 
as shown in Figure 4. The choice of s is not unique; we can 
reproduce visual light curve by adopting another vale of 
s. Here, we choose s to reproduce not only V band mag- 
nitudes but also UV 1455 A continuum fluxes (Cassatella 
et al. 2004). This is a strong constraint for a choice of 
s, and thus, we hardly find another s that reproduce both 
visual and UV light curves. 

Figure 4 depicts our modeled light curve that reproduces 
well both the early optical and UV 1455 A continuum 
light curves. The observed UV flux is small even in the 
super-Eddington phase in which the photospheric lumi- 
nosity is several times larger than that in the later phase. 
This means that the photospheric temperature is as low as 
logT < 4.0. In our model, the temperature is logT = 3.93 
at the optical peak and lower than 4.0 for 8 days after the 
peak, gradually increasing with time. Such a behavior is 
consistent with B — V evolution reported by Chochol et 
al. (1993), in which B — V is larger than 0.3 for the first 
ten days from JD 2,448,677 and gradually decreases with 
time. 

In the later phase, our modeled visual magnitude decays 
too quickly and is not compatible with the observed data. 
Hachisu & Kato (2005) concluded that this excess comes 
from free-free emission from optically thin plasma outside 
the photosphere. They reproduced well the optical light 
curve in the late phase by free-free emission as shown by 



Super-Eddington Luminosity Model 
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the dash-dotted line in Figure 4. 

We see that the peak luminosity exceeds the Eddington 
limit by 1.7 mag, and the super-Eddington phase lasts 12 
days after its peak. 

The distance to the star is obtained from the comparison 
between observed and calculated UV fluxes, that is, 1.83 
kpc with A\ = 8.3E(B-V) = 2.65 for A=1455 A (Seaton 
1979). From the comparison of optical peaks, the distance 
is also obtained to be 1.83 kpc with Ay=0.99 (Chochol 
et al. 1997). This value is consistent with the distance 
discussed by Chochol et al. (1997) that ranges from 1.3 
to 3.5 kpc with a most probable value of 1.8 kpc (see also 
Rosino et al. 1996). 

4. DISCUSSION 

Shaviv (2001b) found two types of radiation- 
hydrodynamic instabilities in plane parallel envelopes. 
The first one takes place when (3 decreases from 1.0 (before 
ignition) to ~ 0.5 and the second one occurs when (3 de- 
creases to ~ 0.1. Here (3 is the gas pressure divided by the 
total pressure. When the luminosity increases to a certain 
value the envelope structure changes to a porous one in a 
dynamical time scale. Radiation selectively goes through 
relatively low-density regions of a porous envelope. Sha- 
viv (1998) estimated effective opacities in inhomogeneous 
atmospheres and showed that they always less than the 
original opacity for electron scattering, but can be greater 
than the original one in some cases of Kramer's opacity. 
Ruszkowski & Begelman (2005) have calculated radiative 
transfer in slab- like-porous atmospheres, and found that 
the diffusive luminosity is about 5-10 times greater than 
the classical Eddington luminosity when the density ratio 
of porous structures is higher than 100. 

In nova envelopes, we do not know either how clumpy 
structures develop to reduce the effective opacity or how 
long such porous density structures last. The exact value 
of the opacity reduction factor s is uncertain until time- 
dependent non-linear calculations for expanding nova en- 
velope will clarify the typical size and the density contrast 
in clumpy structures. Therefore, in the present paper, we 
have simply assumed that s is a function of temperature 
and made it to satisfy the condition that s is larger than 
unity deep inside the envelope and approaches unity near 
the photosphere. 

The anonymous referee has suggested that s may be a 
function of "optical width" over a considering local layer 
rather than a function of temperature. Here, the "optical 
width" means the optical length for photons to across the 
local clumpy layer in the radial direction. If this "optical 



width" is smaller than unity or smaller than some critical 
value, the porous structure hardly develops and then we 
have s = 1. In the opposite case, the porous structure 
develops to reduce the effective opacity and then we have 
s much larger than unity. The "optical width" description 
may be a better expression for the opacity reduction fac- 
tor s, because the relation between the opacity reduction 
factor and porous structure is clearer. 

We have estimated the "optical width" (St) of a local 
layer using the solution at the optical peak in Figure 4: it 
is St ~ 3 near the photosphere, 19 at logT = 4.76, 580 
at logT = 5.56, 2.8 x 10 4 at logT = 6.36, and 2 x 10 7 at 
logT = 8.03, i.e., the nuclear burning region. Here, we 
assume that the "geometrical width" of the local layer is 
equal to the pressure scale hight, r/(d\nP/dlnr). This 
"optical width" decreases quickly outward and reach the 
order of unity in the surface region, i.e., the "optical 
width" is large at high temperature regions and small in 
low temperature regions. Therefore, we regard that our 
assumption of s qualitatively represents the dependence 
of the opacity reduction factor on the "optical width" of a 
local layer. 

In our computational method, this "optical width" is 
calculated only after a solution is obtained after many 
iterations to adjust boundary conditions. The feedback 
from the "optical width" ' requires further huge iterations. 
Therefore, in the present paper, we assume a simple form 
of s. 

The wind acceleration is closely related to spatial change 
of the effective opacity. In a case of varying s, for exam- 
ple, when it is a function of temperature, s determines the 
wind acceleration. If we assume the other spatial form of 
s, the acceleration is possibly very different. In our case in 
Figure 4, s is a monotonic function, and then the reduced 
opacity still has a strong peak at logT ~ 5.2 although 
the peak value is smaller by a factor of s than that of the 
OPAL peak. The resultant velocity profile is essentially 
the same as those in Figure 1; the wind is accelerated at 
the shoulder of the OPAL peak. 

We thank A. Cassatella for providing us with their ma- 
chine readable UV 1455 A data of V1974 Cyg and also 
AAVSO for the visual data of V1974 Cyg. We thank an 
anonymous referee for useful and valuable comments that 
improved the manuscript. This research was supported in 
part by the Grant-in- Aid for Scientific Research (16540211, 
16540219) of the Japan Society for the Promotion of Sci- 
ence. 
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log R (cm) 

Fig. 1. — Velocity (dotted), temperature (T: solid), density (p: solid), local Eddington luminosity (I/Edd : dashed), and diffusive luminosity 
(L r : solid) for three envelope solutions on a I.OMq WD. The critical points of steady state winds are denoted by crosses, circles, and triangles 
for three solutions with s=X, 3, and 10, respectively. The right edge of each line corresponds to the photosphere. 
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Fig. 2. — Photospheric velocity (u p h); photospheric luminosity (Z/ p h)i and wind mass loss rate (dM/dt) for three evolutionary sequences, 
s=l (solid), 3 (dashed), and 10 (dash-dotted) against the photospheric temperature (T ph ). The value of s is attached to each curve. Novae 
evolve from low to high photospheric temperatures. 
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Fig. 3. — Theoretical light curves of a 1.0 Mq WD for reduced opacities of s = 1, 3, and 10. The value of s is attached to each curve. The 
maximum visual magnitude of each light curve is indicated by an arrow. 
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JD 2,448,000+ (day) 

Fig. 4. — Theoretical light curves of our 1.05 Mq WD model for a variable s model: visual (upper solid line) and UV 1455 A continuum 
(lower solid line). Observational visual magnitudes of V1974 Cygni (dot) are taken from AAVSO and observational UV (circle) are from 
Cassatella et al. (2004). Dashed: light curve of 1.05 Mq WD for s = 1. Dash-dotted: visual flux from free-free emission model by Hachisu 
& Kato (2005). Dotted: s is shown in a linear scale between 1.0 (at the bottom) and 5.5 (at the peak). 



